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Abstract 



The pair correlation functions for a bilayer system of classical electrons 
has been calculated through a two-component HNC formalism. The results 
show changes in the structure of the correlation functions as the interlayer 
distance is increased from zero to about twice the Wigner-Seitz radius. These 
changes are indicative of structural phase transitions which are expected to 
take place in the crystallized phase and whose short range orders the liquid 
phases emulate. 

Ms. No PACS numbers: 73.20.Dx, 73.40. Gk. 
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I. INTRODUCTION 



Electronic bilayer systems (2D electron liquids in a neutralizing background, separated 
by a distance d, each having density n = 1/na 2 , a being the Wigner-Seitz radius) have 
attracted a great deal of attention, both in condensed matter physics and in relation 
to the Penning trapped ion experiments It has been realized that correlations 

play a more important role in such systems than in the 2D layer or in 3D bulk systems. 

Although attempts have been made to estimate the interlayer correlation functions ]3],|6],[7] 
no reliable calculation on them has so far been available. In this paper we report on the 
results of HNC calculations of the correlation functions of a classical bilayer for a wide range 
of r = e 2 /akT and d/a values. While in most experimental situations strong magnetic 
fields are present, and in the semiconductor bilayers the electronic liquid is degenerate, the 
classical model is expected to be quite reliable in the high T or high r s domain where the 
electrons are quasi-localized. 



II. CALCULATION 

The calculation is based on mapping the bilayer into a 2D two-component single layer 
with an interaction matrix || 

Mr) = Mr) = ~ (1) 
r 

Mr) = TO (2) 

The mapping makes it possible to apply the two-component HNC formalism, which is well 
known ||. The resulting system of equations is 

9ij(r) = hij(r) + 1 (3) 
hij(r) = exp(Hij - </>ij(r)) - 1 (4) 
Hij(r) = hij(r) -Cij(r) (5) 

M*) - |f (6) 



. / 1,\ c u (k) + A 12 (k)h 12 (k) 

hll{k) = l^A~(k) (7) 

Aij(k) = n Cij {k) (8) 

B(k) = [l-A ll (k)} 2 -[A 12 (k)] 2 . (9) 

The system of equations is solved by using Lado's method [[J with the following parameters: 
maximal radial distance R ~ 29, a ~ 2.8, iV = 500. Using as convergence monitor the 
maximum difference between the input array elements and the corresponding elements after 
one HNC iteration, it was found that the difference for the intralayer correlation function 
always exceeds the difference for the interlayer function. Therefore, only the former was 
monitored, and the iterative solution was obtained for a maximum difference less than 
e = 5x 10" 5 . 



III. RESULTS 

Our results for gu(r) and g\ 2 {r) are summarized in Figs. 1,2 and 3,4,5. For low T 
values the most interesting effect is the generation of a partially filled correlation hole in 
layer 2 above a particle in layer 1. The correlation hole empties both with increasing T and 
decreasing d (see Figs. 1,2 and 7). 

Qualitatively new effects appear for high T values, where the quasilocalization of the 
particles becomes important. It has been proposed [§,[5| that in the crystalline phase (for a 
monolayer T cr i t = 137; for a bilayer the T cr i t (d) has not been established, but probably is of 
the same order of magnitude; see comment below) a number of structural phase transitions 
take place as d/a changed. Here in the liquid phase one sees a rather dramatic manifestation 
of the change of short-layer order which reflects the underlying structural phase transitions, 
as d/a is varied from to 2. 

At d = (when the two layers collapse into one single 2D layer) the two correlation 
functions gn(r) and gi 2 (r) are identical, since the two "species" are now one and the same 
and the correlation function is that of a 2D liquid with double density (2n). In the underlying 
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crystalline phase "species 1" and "species 2" particles are randomly distributed on the 
vertices of a triangular lattice with lattice constant a^ = 1.35a. As d/a is increased gn(r) 
and gu{r) remain similar, but develop different amplitudes, gn(r) < gui^)- This reflects 
the change of the first and second coordination numbers q$ , and q$ , qf^ (estimated 
here from the peak values of rgn(r) and rgi2(r)) and a tendency of the two species to occupy 
alternating rows of the triangular lattice which leads to a centered rectangular lattice (Fig. 9). 
This transition seems to be completed around d/a = 0.2 where the position of the first peak 
of gn (r) and that of guir) begin to separate (Fig. 3), marking the gradual transformation 
of the rectangular structure into a centered square lattice (Fig. 9). At d/a = 0.65 the first 
and the second peaks of gn(r) (R^ and R±x) merge (Figs. 4 and 6) which we can interpret 
as the formation of the latter structure. During these processes the q[^ /qn ratio changes 
from 1 (at d — 0) to its maximum 1.66 (at d/a = 0.4) (Fig. 7), qualitatively corroborating 
the above picture. The abrupt formation of a distant second shell at d/a = 0.65 is indicative 
of an underlying (equilateral) square or rhombic structure. Indeed, for d/a > 0.65 the 
transformation into a centered rhombic lattice (Fig. 9) can be seen by the shrinking of 
(Fig. 6). By d/a = 1.5 the rhombic structure transforms into a centered triangular lattice 
(Fig. 9) where the particles of the second layer sit over the holes created in the first layer. 
This is indicated by the reversal of the decrease of (as particle 2 shifts from the center 
of the rhombus into the center of the triangle, and the settling of Ryy , R\i and R^ near 
their expected values a^, y/3a& and , with a/\ = 1.90a, respectively. Although the 

staggered triangular is the final equilibrium structure for d/a > 1.5, further increase of d 
weakens the correlations between the layers to the point that by d/a = 3 the interlayer 
correlations practically vanish (Fig. 8). 

The d/a values of the structural phase transitions predicted from the comparison of the 
respective Madelung energies show reasonable agreement with our results. However, it 
should be emphasized, that in the considered T-range system is liquid and the various lattice 
sites are not actual positions of the electrons, but rather positions of local energy minima, 
and indicative of the positions of the succesive high probability "shells" in the liquid phase. 
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Finally, one may wonder whether the bilayer system crystalizes at a lower T crit value 
than a single 2D layer, as conjectured by Ref. 3. Although no definitive answer can be given 
on the basis of the calculations performed here, an inspection of Figs. 1 and 2 shows that 
there is certainly no tendency for the short-range order to set on a lower T-values than in 
the case of a 2D layer. 

(This work has been partially supported by NSF Grant PHY-91 15714.) 
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FIGURES 

Fig.l Intralayer correlation function for different T values at d/a = 1. 
Fig. 2 Interlayer correlation function for different T values at d/a = 1. 
Fig. 3 Intralayer and interlayer correlation functions. 
Fig. 4 Intralayer and interlayer correlation functions. 
Fig. 5 Intralayer and interlayer correlation functions. 

Fig. 6 Positions of the first and second peaks of the intralayer and interlayer correlation 
functions as indicators of the nearest neighbour positions 

Fig. 7 Normalized values of the first and second peak amplitudes of the interlayer and 
intralayer correlation functions as indicators of the first and second coordination numbers, 
q(l and gg. 

Fig. 8 The value of the interlayer correlation function at r = for different T values, 
indicating the degree of correlation between the layers. 

Fig. 9 The four principal lattice structures. 
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http://arXiv.org/ps/cond-mat/9701195vl 



